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Nuclear problems in astrophysical <7-plasmas and environments 
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Experimental measurements in terrestrial laboratory, space and astrophysical observations of variation and 
fluctuation of nuclear decay constants, measurements of large enhancements in fusion reaction rate of deuterons 
implanted in metals and electron capture by nuclei in solar core indicate that these processes depend on the 
environment where take place and possibly also on the fluctuation of some extensive parameters and eventually 
on stellar energy production. Electron screening is the first important environment effect. We need to develop 
a treatment beyond the Debye-Hiickel screening approach, commonly adopted within global thermodynamic 
equilibrium. Advances in the description of these processes can be obtained by means of g-thermostatistics 
and/or superstatistics for metastable states. This implies to handle without ambiguities the case ^ < 1. 
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1. INTRODUCTION 

Many physical laws are established in the frame of the hy- 
pothesis that each single event of a component of an observed 
system is independent on any other event. For instance, the 
nuclear decay law is based on the hypothesis that each nu- 
cleus decays without being influenced by the other nuclei and 
that all nuclei have the same probability of decaying. 
If P{dt) = Xdt is the probability of decay of a nucleus in the 
time interval dt, the probability of surviving after n intervals 
dt is given by the power law 



(1 -P{dt))" = {l-Xdlf = (i 



(1.1) 



By making n going to infinity the exponential law exp(— A,f ) 
is obtained |T, 7, 3,7]. 

Many physical phenomena are studied by assuming that all 
typical events are independent and non correlated as, for in- 
stance, the nuclear reaction events that occur to form nuclear 
fusion rates in stellar cores. 

Therefore, one uses the Maxwell-Boltzmann (MB) distribu- 
tion to describe both the ionic and electronic component. At 
most, one assumes that all effects from particle correlations 
and/or non linear effects can be neglected, being their contri- 
bution small. 

Also Debye-Hiickel (DH) approach to electron screening, de- 
veloped to take into account the electron and ion influence 
over the pure Coulomb potential of a given ion charge, is 
based over additive, linear effects and standard exponential 
distributions nil Sil. 

Very recently, in many different experiments, influence and 
effects of environment have been observed in addition to the 
known effects over fusion thermonuclear reactions in stellar 
cores. One can assume that these effects are due to cor- 
relations, microfield distributions, electron screening, ran- 
dom forces, fluctuations of extensive quantities, among others 
causes. Here we give a list of the recent observations: 

1 . Measurements of electron capture decay rates and con- 
stants of ^Be in host materials iBIHElIilB, of a- 
radioactive uranium exposed to glow-discharge plasma 



fTT] and of other heavy nuclei implanted in different 
environments [15, 16, 17] have recently shown fluctua- 
tions and variations of previous standard results. 

2. Rates of deuteron-deuteron fusion in metal matrices 
have shown a great increase respect to the standard lab- 
oratory measurements and have shown a non negligible 
influence of the conducting electrons LIS^ ,19> ,20^ 

3. Variations observed in electron capture by ^''K could 
be responsible of observed discrepancies in the density 
ratio U / Pb in Oklo site phenomenon ll23ll . 

4. Unexpected and unexplained fluctuations in the decay 
rates of ^^Si and ^^^Ra have been reported and evidence 
of correlations between nuclear decay rates and Earth- 
Sun distance has been found (Jenkins-Fishbach effect 

illlllli)-^ 

5. Several careful experiments designed to study the de- 
cays of long lived radioactive isotopes have reported 
observations of small periodic annual variations mod- 
ulating the well-known exponential decay curve. If 
Jenkins-Fishbach effect is correct we could have to con- 
sider profound consequences for many areas of physics 
and engineering. Discrepancies in half-life determina- 
tions of many nuclides reported in literature have been 
attributed to variations of solar activity during experi- 
mental works or to seasonal variations in fundamental 
constant. However, while in a satellite experiment the 
correlation between variation of rate and Earth-Sun dis- 
tance has not been observed ll28ll . other authors suggest 
that neutrinos from the flare of 13 December 2006 have 
induced a change in the decay rate of ^^Mn. 



Very recently, after new measurements, the variation has been reduced and 
an experimental group suggests that non correlation exists with the Earth- 
Sun distance, although modulation and variation still exist 12711 
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6. Precise half-life of unstable isotopes, in particular 
of '^C, are considered in doubt, since the Jenkins- 
Fishbach effect implies that the decay rate of an iso- 
tope is not a constant of Nature, in contrast with the 
findings of Rutherford, Chadwick and Ellis in 1930 that 
the rate of transformation of an element is constant un- 
der all conditions. A reasonable interpretation of the so 
called wiggles in ^'*C decay rate has been suggested by 
Sanders: wiggles indicate that the solar fusion furnace 
is pulsating in some way like the Cepheyd furnace i29ll . 

7. Recent work by Sturrock on power spectral analyses 
of radiochemical solar neutrino data and solar irra- 
diance have revealed modulations attributable to non 
spherically-symmetric and variable nuclear burning in 
the solar core [30]. 

8. New determinations of solar CNO content based on 3D 
hydrodynamic model and a new treatment of departure 
from Local Thermodynamical Equilibrium (non-LTE) 
have mined agreement between helioseismology and 
observations IBlll . 

9. Local non-Gaussianity of the temperature anisotropy 
of the cosmic microwave background and other non- 
Gaussianities revealed by satellite experiments have 
been discussed recently in many papers to which we 
send for details, since these problems are not object of 
our interest here. 

Of course aU the above effects are not concerning with the 
normal evolution of a star, rather represent fluctuations around 
average values, the time averages giving rise to constant val- 
ues. The states at which the nuclear processes take place are 
not standard global thermodynamic equilibrium states. They 
can be seen as metastable states with a given life time (eventu- 
ally a long-lifetime) and therefore the use of ^-thermostatistics 
1I32I1 and/or superstatistics ll33i [34ll can be the appropriate ap- 
proach to be applied to understand the above effects ll35ll . 
In most of all these processes, the first and most important 
environment effect that operates is the electron screening. 
Usually it is described by DH potential: Coulomb potential 
which originates from a charge in an astrophysical plasma is 
screened by the other charges. Decrease of potential is faster 
than pure Coulomb one because of the exponential factor in 
DH approximation. 

To have a more complete description one must go beyond DH 
potential approach; the electron cloud that screens can be, for 
instance, more concentrated around the ion that originates the 
Coulomb potential and its momentum distribution can differ 
from the MaxwelHan one lisi I2lll36ll . 

Nuclear fusions in stellar cores are influenced by electron 
screening. Among many nuclear astrophysical processes 
that are influenced by the surroundings, the electron capture 
by nuclei, like for instance e^+^Be^^Li+v in the Sun or 
e^ +p^n+v in high density stars, plays a special role because 
the electron component of the plasma is responsible of the 
screening and the same electrons are themselves responsible 
of the screened physical process that is the electron capture 
by nuclei from bound and continuum states caused by weak 



nuclear interaction. The first action is due and characterised 
by the charge distribution in the coordinate space, with this 
distribution we can build the screening potential. The second 
action is the weak nuclear process of capture with emission of 
a neutrino. Its rate is characterised by the momentum distri- 
bution of electrons in the momentum space. 
Why we are looking for a screening potential different from 
the widely used DH one? 

If it is required that the distribution in coordinate space be 
more concentrated around the ion than in DH case, then this 
can be realized by a ^-distribution with q < 1 (depleted tail). 
For instance, we can show in the case of electron capture by 
^Be, that the electron distribution tail should be depleted and 
this fact agrees with the requirement that to obtain an increase 
of the EC rate (to obtain a decrease of the ^B neutrino flux). 
The number of particles inside the DH sphere is, for many 
astrophysical systems like stellar core, very low. Therefore 
standard MB statistics does not apply. The fluctuation of den- 
sity is large and due to the connection of density with tem- 
perature through the equation of state also the fluctuation of 
temperature is present. We may use superstatistics to justify 
our choice of non extensive statistics and its distribution. 
Coulomb cutoff in the spatial space was already used. This 
is the simplest way to screen the Coulomb field: to cut it at a 
certain r — R so that the potential is zero outside R. No other 
justifications have been advanced for this approach. 
A result analogous to this can be obtained considering den- 
sity and temperature fluctuation in superstatistics or non lin- 
ear Poisson equation, in such a way that use of the non exten- 
sive Tsallis statistics with ^ < 1 is justified. Of course q < I 
means that the electron spatial distribution is more concen- 
trated near the positive ion than the distribution responsible of 
DH screening which is the Coulomb potential times an expo- 
nential factor 

Moreover, we have two coordinate spaces to consider: spatial 
coordinate and momentum coordinate. If <7 is a universal pa- 
rameter for the observed system the same q should be used for 
the spatial and for the momentum (energy) coordinates sys- 
tem. We assume that the mechanism that motivates a ^ < 1 
distribution, due to the density fluctuation, is also responsi- 
ble, through the relation density-temperature-momentum (en- 
ergy) of the momentum q < I distribution. Therefore also 
the momentum component of the distribution is a non exten- 
sive Tsallis distribution with q < I. On the other hand, if the 
same q occurs for spatial coordinate and for momentum co- 
ordinate distributions, this could imply some problems with 
uncertainty principle with ^ < 1 systems |37tl- 
Other requirements and observations impose to use a <7 < 1 
distribution in the momentum space. 

For instance, in the study of the spectral lines of the Sun a 
cutoff is required. As shown, for example in |'38','39'] the elec- 
tron component of a central core stellar plasma is made by 
correlated electrons where fluctuations of several types occur. 
The equation of state contains many non standard terms and 
we can realize that in such a case statistical distribution must 
deviate from MB one. 

Although very recent experiments at GSI concern with oscil- 
lation of decay constants (known as GSI anomaly) we do not 
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discuss here this case between the environment effect could be 
absent due to the specific experimental conditions [I40[l4lll42ll . 

In this work we want to limit ourselves to discuss only two 
nuclear problems of astrophysical interest where we use the 
Modified Debye-Hiickel (MDH) approach. Before this dis- 
cussion we describe how one can go beyond the DH approach 
and derive a screening potential in an astrophysical plasma 
using a non extensive q < \ electron distribution (Section 2). 
As an application of the MDH potential we study the d-d fu- 
sion rate in metal matrices. Also in this case we take q <\ 
and show that a good agreement with the experimental results 
can be obtained (Section 3). The electron capture decay by 
^Be is discussed, the variation of the rate evaluated and few 
implications reported about the neutrino fluxes (Section 4). 
Then we discuss the meaning, on a microscopic level, of tak- 
ing ^ < 1 for particle distributions. Formation of a distribution 
with ^ < 1 is discussed and conclusions are reported (Section 
5). 



2. MODIFIED DEBYE-HUCKEL SCREENING 

Thermal effects and screening phenomena in plasma envi- 
ronment are strictly connected. Assuming a MB thermal dis- 
tribution, the DH ion-screening potential 



y..M = ^exp(^-— j , 
is obtained through the Poisson equation, where 



R 



kT 



DH 



4-nnZpe^ 



(2.2) 



(2.3) 



is the DH radius, but non-local thermodynamic effects can de- 
viate the high energy tail of the velocity distribution from the 
exponential MB form. 

Recently f^, it has been proposed an approach based on the 
application of super-statistic and/or a ^-version of the Poisson 
equation to the DH screening model by assuming that non- 
linear effects produce fluctuations on the inverse of the DH 
radius 1/Rdh, with a Gamma-function probability distribu- 
tion 

Mr, X, Xo) = ^li!jMIlxThr^ e-^A,irM) , (2.4) 

where fg{r,X,Xo) represents the probability density to ob- 
serve a certain value X spreads around a central value A,o = (X), 
with 



(X)^ / f,ir,X,Xo)XdXi 



1 



Rdh 



(2.5) 



To obtain from fg{r,X,Xo) an electron depleted tail distri- 
bution with a cutoff, we limit the entropic index q into the 
< ^ < 1 interval. Furthermore, we assume 



Xo) 



1 



(2.6) 



where g{q) is a generic entropic index function that satisfies 
the condition g(l) = 1 (in L8i,,36J the choice g{q) — \/{2~ q) 
is adopted). 

The point charge potential Vq{r) can be identified with the 
functional 



J,(r,Xo) = j f,{r,X,Xo)e-^' 



dX. 



(2.7) 



through the relation 

rVci{r)=ZiZ2e^ 1^-^^ Jc{r,Xo) ■ (2-8) 

The charged particles distribution p^(r) and the point charge 
potential Vq{r) oc p^(r) can be derived from Eq. ( 12.8b . devel- 
oping the functional fq{r, Xq) with the previous assumption. 
A cutted form for the potential is obtained as 



exp^ 



where 



and 



-1 



exp^(x) = [l + (l-^)x]y('-'') 



J+ > 



(2.9) 



(2.10) 



(2.11) 



is the ^'-exponential function with — x for x > and 
= for X < 0. In this way, Yqir) vanishes for r > 

The main contribution to the charged particles fusion cross 
section is given by the screening barrier penetration factor 
P{E). In the standard DH potential case, the simplified ex- 
pression 



(2.12) 



can be obtained iS{E) is the astrophysics factor), that differs 
from the bare nuclei cross section Obare — S{E) /E only for the 
penetration factor 



P{E) = exp -71 



Eg 



E + Udh 



where 



D 



DH 



(2.13) 



(2.14) 



with D = Zi Z2 and Eg is the Gamow energy. 
Differently, in the case of the MDH potential Vq{r), the pene- 
tration factor Pq{E) is given by the expression: 



Pq{E) = exp j I ^2ijc^ (Vqir) -E)dr\ , (2.15) 
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where the classical turning point rg has to be determined 
through the equation Vq{rg) — E. 

Equation ( 12.151 ) can be analytically solved in the ^ = case, 
obtaining 



D 



E + g{0)D ' 



and 



Eg 



E+g{0)D_ 



(2.16) 



(2.17) 



The bare nuclei cross section can be corrected, to account for 
the MDH screening, multiplying Ob^reiE) by the factor 



Pg{E) exp 71 




(2.18) 



In conclusion the penetrating charge gains the energy Udh in 
the DH approach and the energy g{Q) D in the MDH approach 
when ^ = 0. 

Let us remark that the function g{q) is not arbitrary but is 
related to the charge distribution p^. In fact, it is easy to show 
that for ^ < 1 



(2.19) 



so that the screening charge is distributed from r = to 
'"cut = ~ l)^<:j]- Therefore if we have information from 
experiments or from models on the charge distribution we 
can fix the function g{q) and its value. For instance, in 1I22I1 
the distribution around deuterons is evaluated by means of a 
Thomas-Fermi model. 



r'/^, but the mean number of quasi-free particles in the De- 
bye sphere results much smaller than one. Then, the picture 
of the Debye screening, which should be a cooperative effect 
with many participating particles, seems not to be applicable. 
The thermal motion of the target atoms is another mech- 
anism capable of increasing the reaction rate; however, 
Maxwellian momentum distribution at the experimental tem- 
peratures gives negligible effects. The relationship between 
energy and momentum of quasi particles can be broaden by 
many-body collisions, then a long tail, non-Maxwellian mo- 
mentum distribution can emerge from a MB energy distribu- 
tion. Fusion processes select high-momentum particles that 
are able to penetrate the Coulomb barrier and are, therefore, 
extremely sensitive probes of the distribution tail ll43l l44ll . 
This quantum dispersion effect has been already introduced 
as a possible explanation of the reaction rate enhancement. 
The screening plasma particle effect on the reaction rate has 
been evaluated adopting for the first time, the MDH potential 
proposed in fst]. 

For instance we compare the MDH astrophysical factor for 
d-d reaction with the experimental data reported in ifTol I20I1 . 
adopting the choice g{E) = 3 — 2q. The results is shown in 




3. ANOMALOUS ENHANCEMENT IN LOW ENERGY D-D 
FUSION RATE 

In recent years a number of different experiments, with tar- 
get adsorbed in a metallic matrix, have evidenced a strong en- 
hancement in the fusion reaction rate at few keV |18, 19, 20, 
lUll^l. For instance, the d(d, t)p reaction has been widely in- 
vestigated and ^•^Li(d, a)^'^He has been studied with similar 
results. 

A less strong enhancement has been observed in gas target ex- 
periments, that can be easily explained by the standard elec- 
tron screening, with a potential 11^ of the same order of the 
adiabatic limit U^d = 28 eV. 

However in deuterated metal target experiments, a potential 
Ue of hundreds of eV, ten times greater than the limit Uad, is 
needed to reproduce the results. 

A p ossible explanation has been proposed (see ifTsl [T9I, l20l 
ill Hi) based on a simplified model of the classical quasi- 
free electron. It predicts an electron screening distance of the 
order of the Debye length. 

This approach reproduces both the correct size of the screen- 
ing potential Ue and its dependence on the temperature: 11^ °^ 



Ql , , . , , 1 

0.5 1 2 5 10 20 



FIG. 1: Astrophysical Factor experimental points from ref. fT^l20ll . Bare 
nuclei curve correspond to Sbaie = 43 + 0.54£'cm keV b, while the screened 
curves are 5 = /, ■ 5ba,c(£cm), = (I/Rdh), giq) = 3 -2q and q = 0. Ion 
thermal motion is neglected: E^m = E/2. 



Fig. [T]for the entropic index q — Q. 

One can observe as a screening potential Ug three times 
lower than in the standard DH case (Ug ~ lOOeV instead 
of Udh ^ 300ey), is required to reproduce the experimen- 
tal data. An electrostatic screening potential of this order of 
magnitude has been obtained, for instance, by Saltzmann and 
Hass 1 22] through a Thomas-Fermi model of the electron gas 
in a deuterated-copper target (they obtained a screening po- 
tential of 163 eV, instead of the 470 eV needed to reproduce 
the experimental results). 

Now we are treating g{q) and q as free parameters, but, in 
principle, a link can be establish between the inverse DH ra- 
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dius, the temperature fluctuations and the q'-index: 



^(]JR^^^(kT) 

I/Rdh kT 



q. 



(3.1) 



By this way the modified DH potential can be obtained start- 
ing from the environment condition. 

We conclude: the reaction cross section enhancement, due 
to plasma screening effects, has been evaluated adopting the 
MDH potential. As shown in ifisi [T9I1 . accounting the ma- 
trix metal valence electrons, a screening potential U w lOOeV 
can be obtained, greater than the adiabatic limit U^d but even 
too low to reproduce the experimental results. We showed as 
this discrepancy can be removed adopting the MDH potential 
( 12. 9t with a proper choice for the function g{q) and the en- 
tropic index q, as can be seen in Fig. [1] In principle the value 
of the entropic index q can be derived from the plasma equa- 
tion of state, then the MDH screening contribution to the rate 
enhancement can be evaluated exactly. We will investigate 
this last point in a future work. 



4. ELECTRON CAPTURE BY ^BE 

Electron capture (EC) is a nuclear process where an inter- 
nal proton of a nucleus is converted to a neutron by means of 
weak interaction: e" +(Z+ 1 , A)^(Z, A)+v. 
In laboratory, this reaction proceeds by capture of a K-shell 
electron since these electrons have the greatest probability to 
move in the nuclear region. In the stellar core, atoms are al- 
most completely ionised so that EC proceeds with continuum 
electrons and the rate depends on the electron density over the 
nucleus. In plasmas of very high density, stable nuclei are 
forced to capture free energetic electrons and hydrogen gas 
change in a neutron gas. 

In stars like Sun, the capture e^+^Be^^Li+v is one of the 
most important p-decay reaction. In a star in nuclear equi- 
librium, p-decay plays a special role because the inverse re- 
action requires the absorption of a neutrino. The neutrino 
escapes from the local environment with an average energy 
E^, = 0.814Mey. Therefore, p-decays cannot participate to a 
true equilibrium; they are sufficiently slow so that an equilib- 
rium can be established in a time short compared to the time 
required for a significant change of the element average com- 
position. However, environments must be such to satisfy this 
time condition and this can be so if only frequent and pure 
Coulomb interactions are active. 

The principle of detailed balance ensures that in global ther- 
modynamic equilibrium loss of electrons captured by nuclei 
in the medium-high energy tail is counterbalanced by produc- 
tion of such electrons due to inverse reaction. However, neu- 
trinos, once produced, are travelling away with their kinetic 
energy and the intermediate energy region of electron distri- 
bution results depleted compared to Maxwellian distribution. 
In fact, the main contribution to the EC rate is from the region 

at 2 3 time kT = 1.21 keV. Electron component is a non- 

LTE state. This metastable state can be well described by a 
^-distribution function with ^ < 1 . 



It is known that in nuclear continuum EC the rate includes 
the Fermi factor, i.e. the electronic density at the limit of r ^ 
in a pure Coulomb potential. The pure nuclear cross section 
is corrected by a factor because of the Coulomb interaction 
between the captured electron and the nucleus. Usually, in a 
plasma environment, DH potential takes the place of Coulomb 
potential. In this case the electron density at the nucleus is 
known only numerically from the solution of the appropriate 
Schroedinger equation. 

In an astrophysical plasma environment, the momentum dis- 
tribution of screening electrons can differ significatively from 
MB distribution. Here we use the MDH potential described 
previously. 

Of course, also in this case we need to evaluate the electron 
density at r = by solving the appropriate Schrodinger prob- 
lem for electrons in the continuum. 

Whereas the solution with a Hulthen potential can be given 
in a close analytic form, the MDH potential (as well as the 
standard DH potential) admits only numerical solutions for 
the electron density at r = 0. Alternatively, one can use the 
Hulthen potential that fits quite well the MDH potential in the 
small r region (near the nucleus), but contains an infinite tail 
instead of a cut-off at an appropriate value of r = rcut, as for 
the MDH potential with q<l. 

In the following, we evaluate the rate for the free electron cap- 
ture by a (A,Z)-nucleus, given by the integral, in the three 
dimensional space of velocities, of electron capture cross- 
section CTe times the electron velocity v, the normalised proba- 
bility density {F^ , Fpjj , F^^oi F^)^ and the normalised velocity 
distribution for electrons given by the function 



'^'ilKkT) 



3/2 



exp„ 



(4.1) 



where q = I for C, DH and H (with /^^^ (v) = /^g (v) the nor- 
malized MB distribution).^ We define the pure Coulomb nu- 
clear electron capture rate, averaged over a MB distribution, 

as 

^jr) = I {a,v)F^{E)f^^iv)4%v^dv, (4.2) 



where 



a. = 



c 



(wo+wy 



(4.3) 



is the nuclear electron capture cross section |36|], with G the 
Fermi constant. Wo the nuclear energy release for one electron 
with total energy W,x^Cl{l)^ + Cl{af the well-known 



- The normalised probability density that an electron of the continuum spec- 
trum with velocity v and travelling in a screening potential (Vc, Vdr, Vh or 
Vq) be at the nucleus with coordinate r = 

^ For details on this approach, expressions used and complete numerical re- 
sults see reference I3q|. 
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reduced nuclear matrix element. 

The Fermi factor for Coulomb potential, given by 



271 ri 



1 



,-271 r| 



(4.4) 



with r| = 4/ (flo p), where is the Bohr radius and p — nig v 
is the electron momentum, follows from the definition 



FciE) 



lim 



¥c('-) 



pr 



(4.5) 



where \\t^ (r) is the wave function of the Schrodinger equation 
with the Coulomb potential. 

The non extensive rate ^ g, can be obtained by substituting in 
Eq. ( 14.21 ) the factor (E) with the Fermi factor 



Fg{E) = lim 



pr 



(4.6) 



where \\ig{r) can be obtained as a numerical solution of the 
Schroedinger equation with the MDH potential iQ] 



(4.7) 



and^(^) = 1/(2-^). 

Consistently with the derivation of the MDH potential in the 
definition of !^ ^ we must insert in place of the MB distribution 
the normalised non extensive distribution fq{v). The integral 
for the rate l{_q, when q< 1 , is performed over the real interval 
[0, Vcut] with 



2kT 



Vcut — 



(1 -q)mt 



<1 , 



(4.8) 



which defines a cut-off condition in the velocity space. When 
\ the rate %_q reduces to DH rate ^py. 
Although electron density at r = due to Vq{r) is smaller 
than Coulomb density, in the velocity space the probability 
density in the low momentum region is greater than MB be- 
cause the continuum electron distribution fq{v) we use privi- 
leges low momentum electrons. Therefore, screening may be 
important in continuum EC rate. We have calculated devia- 
tions of the rate 1{_^{T) (with X = DH, H and q) respect to 
llciT), SLtkT ^ 1.21 keV (where EC by '^Be takes place). 
For any q < 1, ^q > ^f^. We have verified that deviations de- 
pend very smoothly on kT except for t, < 0.45 ao and depend 
very strongly on q. 

The value of q for EC by ^Be in solar plasma can be de- 
rived from the expression that links q to fluctuation of l/R 
ll36ll . In the solar core, where the average electron density 
is He = 9.1 and the number of particles inside the Debye 
sphere N^^^ is about 4, we can obtain q = 0.86. It is more safe 
to consider a range of values of q between 0.84 and 0.88. At 
kT = 1.21 keV and ^q = 0.45flo the calculated ^tqiT) is esti- 
mated to be about 7 - 10 % larger than standard DH {q=l) es- 
timate that is the 0.69% smaller, at the same conditions, than 
Coulomb rate ^^,{7). Of course, a smaller value of q should 



imply a much greater enhancement of EC rate over DH one. 
Let us consider the ^Be - p fusion. This reaction producing **B 
and, as a consequence, ^B neutrinos, in competition with ^Be 
electron capture. We have verified that the effect of the MDH 
potential over its rate is negligible. In fact, correction to F^ is 
effective only at relative ^Be - p energies lower than 2.4 keV 
where fusion cross section has a negligible value because its 
most effective energy is at l^keV. Therefore, if EC rate of 
^Be increases over its standard evaluation of a given percent- 
age, ^Be increases its destruction while the neutrino flux from 
^Be does not change because the ^Be density decreases. How- 
ever, the ^B flux should diminish of the same percentage. This 
behaviour is in line with what is found in experiments. 



5. PROBLEMS ABOUT DISTRIBUTION WITH g < UTS 
FORMATION AND CONCLUSIONS 

Among the many astrophysical problems and experimen- 
tal observations listed in the Section 1, reporting deviations 
from the standard behaviour detected and measured"* that, to 
our opinion, showing how the environment effects could be 
possibly taken into account by means of g'-thermostatistics 
and/or superstatistics approaches, we have selected two ex- 
amples: the d-d fusion in metal matrix and the solar EC by 
^Be. 

Measured rates of the first subject are strongly enhanced re- 
spect to standard laboratory measures. We have proposed an 
explanation based on the ^ = thermostatistic description of 
the MDH potential used for the d-d fusion reactions. The pa- 
rameter q ^0 means that a great deformation of the standard 
MB momentum distribution function should result evident. 
In the second subject, the electrons that screen the capture 
have a spatial distribution with q — 0.84 — 0.88. Although in 
this case deformation is small, consequences in the production 
of ^B solar neutrinos are of great interest. 
In both cases, the parameter q is lower than 1 and in the sec- 
ond subject the ^-approach is used for the derivation of the 
MDH potential in the coordinate space and for the calculation 
of the rate in the momentum coordinate space and the same 
value of q is assumed. 

We have tried to show how the case ^ < 1 be important 
in many nuclear physics problems in a non extensive q- 
environment. Of course also the case ^ > 1 has important 
applications when, for example, gamma rays produced in the 
plasma comptonise electrons making a fat distribution tail. 
However, to handle with the ^ < 1 case means that we have to 
pose attention to some problems that we list below: 

1 . When using superstatistics to derive ^ < 1 distribution. 
Because the average value depends on a variable, a sec- 
ond average is needed; 

2. The function g(^) of expression l/{2 — q) or 3 <7 — 2 (see 
the derivation of MDH), which goes to 1 for q going 



■ for some of them results should be confirmed. 
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to 1, is arbitrary and we must find a condition for its 
determination. This can be accomplished if we know 
the electron charge distribution around the ion, as we 
have already discussed in Section 2; 

3. When the process requires a double use of non- 
extensivity (we have the system in phase space) in the 
spatial coordinate and in the momentum coordinate, we 
must establish rules in order to know if the same q 
should be adopted or two different values and if they 
are possibly linked in some way, i.e. if an uncertainty 
principle holds; 

4. Particular experimental results are modulated in time; 
in this case one needs a function q{t) rather a parameter 
a- 

We want now to show how ^ < 1 distribution can be formed 
starting from MB distribution. 

The ^ < 1 distribution has a depleted tail with a cut-off and 
an enhanced head. To understand how this distribution can 
be formed, let us start by considering a system of nonin- 
teracting particles distributed along a Maxwellian shape, at 
temperature T . 

We introduce a quantity called cut-off energy £cut, analogous 
to the Fermi energy level of quantum distributions, located 
at Ecut = ^7^/(1 ^?). Depending on the value of q, with 
< q' < 1, the cut-off energy may assume a value in the range 
kT < Ecut < °°- 

The number of particles in the Maxwellian distribution 
with an energy e above the cut-off energy, i.e. within 8cut < 
e < oo, can be easily calculated by 



2 /3 1 
y/K \2 \-q 



(5.1) 



where r(a, x) is the inverse incomplete gamma function. 
We assume that at a certain time a many-body interaction be 
active among the particles. As a consequence of complete oc- 
cupancy principle, all the particles occupying the states above 
8cut displace below leaving all states above empty. Then, par- 
ticles reorganize collectively their distribution at the appropri- 



ate temperature T' <T, where T' 

R 



j^jj^ T spending the work 



L 



ci<l- 



The complete energy balance of this ideal process producing 
a depleted distribution is given by 



where Y(a, x) = r(fl) — r(fl, x) is the incomplete gamma func- 
tion. Since generally 5A^+ /N is negligible, we can write 



-NkT 
2 



where 



' 5-3q 



(5.3) 



3^/K \2' l-q 



2' l-q 



(5.4) 



We remark that L^^^/NkT is always negative because the 
work is done in favor of the environment. It depends on q 
only and shows a maximum at ^ « 0.7. 
We note that MB distribution of non interacting particles in 
thermodynamical equilibrium can be viewed as a generalized 
distribution with q < \ and quasi-particles below £cut dis- 
posed above if the energy level sequence remans unchanged. 
We remember that q < \ case is related to a finite heath bath 
as shown in the past by Plastino and Plastino l l45i l46ll47ll . 
In recent years, many experiments have shown results that de- 
viate from the estimations evaluated on the basis of global 
thermodynamic equilibrium, of independence of single events 
neglecting particle correlations, fluctuations and non linear 
effects. Under these experimental evidences we realize that 
many processes that were understood on the basis of simplifi- 
cations can be better understood if we improve earlier evalua- 
tions by generalising previous approaches. In such situation 
we are convinced that ^-thermostatistics will play a funda- 
mental role. 



<7<1 ' 

(5.2) 
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